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1. Introduction 

The Einstein-Cartan theory of gravity (EC) extends Einstein's theory of general 
relativity (GR) in a natural way by including the spin properties of matter and their 
influence on the geometrical structure of space-time. By removing the symmetry 
requirement on the two lower indices of the connection, Cartan [1] showed that the 
dynamics is no longer entirely determined by the metric; the antisymmetric part of 
the connection called torsion became an independent dynamical variable. Besides 
the energy-momentum of the matter content sourcing curvature, its spin was later 
postulated to be the source of torsion |2J. The EC theory locally satisfies the Poincare 
symmetry [3] accounting for translational degrees of freedom associated with curvature 
and rotational degrees of freedom linked to torsion. 

Weyssenhoff and Raabe initiated a careful study of the behaviour of perfect fluids 
with spin jl] . In order to build cosmological models based on the EC theory, Obukhov 
and Korotky extended their work [5j. They showed, in particular, that by assuming 
the Frenkel conditional the model reduces to the description of an effective fluid in 
GR where the effective stress-energy momentum tensor contains some additional spin 
squared terms. 

The first studies of perturbations of a perfect fluid within GR were carried out by 
Lifshitz [6] in a fixed gauge and reformulated in terms of gauge-invariant variables by 
Bardeen [7] . The dynamics of such a fluid have also been investigated in a more physical 
and transparent gauge- invariant manner by Hawking [8] and extended by Ellis [9]. We 
shall follow the latter approach here and use the 1 + 3 formalism. 

As Puetzfeld points out [10], there are an increasing number of theoretical reasons 
for studying cosmological models based on a non-Riemannian geometry, as some key 
features of the current concordance model such as dark matter, dark energy and in 
particular inflation still need to be explained. The Weyssenhoff fluid, for example, seems 
a promising candidate to describe cosmological inflation in a geometrical manner without 
using scalar fields, which have not yet been observed. This promising behaviour may 
arise from the spin density squared terms contained within the effective stress energy 
momentum tensor derived by Obukhov and Korotky [5], since these spin contributions 
dominate the dynamics at early times. Although the Weyssenhoff fluid is expected to 
leave the late time dynamics unchanged, making it an unsuitable candidate to describe 
dark energy, it may still therefore significantly affect the early time evolution of the 
fluid. 

In this publication, we restricted our study to the formal derivation of the dynamical 
relations for a Weyssenhoff fluid. A detailed study of the large scales dynamics of 
such a fluid in an attempt to get a spin based inflation will be pursued in further 
work. To remain as general as possible we chose not to perform a first- or second- 
order perturbation analysis for a particular class of models. This can easily be done 

| Note that the Frenkel condition arises naturally when performing a rigorous variation of the action. 
It simply means that the spin pseudovector is spacelike in the fluid rest frame. 
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according to the symmetries of the models, and some specific examples will be pursued 
in a later publication. The derivation of the Weyssenhoff fluid dynamics is a prelude to 
the perturbation analysis, which is especially relevant to study the structure formation 
seeded during the inflationnary era. The dynamics of such a fluid in a 1 + 3 covariant 
approach has been studied previously in a cosmological context by Palle [TT]. However, 
the use of effective GR relations in conjunction with EC identities is rather opaque in 
this work, and also certain length scales are excluded from the analysis making a new 
study, which considers all length scales, appropriate. 

In the standard GR theory, the 1 + 3 covariant approach leads to six propagation 
equations and six constraint equations. These give respectively the time and spatial 
covariant derivatives of the set of dynamical variables, which are the energy density p, 
the expansion rate G, the shear density a, the vorticity density u, the 'electric' part 
of the Weyl tensor E and the 'magnetic' part of the Weyl tensor H. The Weyssenhoff 
fluid is described by an effective GR theory, where the additional degrees of freedom 
due to torsion are entirely determined by the spin density S. Therefore, in addition to 
the spin density modifying the dynamical equations for the six standard variables, we 
also expect to find additional dynamical relations. 

In the next section, we briefly outline the EC theory, then give a concise description 
of a Weyssenhoff fluid in Section [31 Section 0] is devoted to the Weyssenhoff fluid 
dynamical analysis using the 1 + 3 formalism outlined in Appendix A. The consistency 
of the particular case with zero vorticity and peculiar acceleration (u = a = 0) is 
established by evolving the constraints in Section [5j The last section draws a comparison 
with Palle's results. In this paper, we use the (+, — , — , — ) signature. To express our 
results in the opposite signature used by Ellis [32] , the correspondence between physical 
variables can be found in [13] and in Appendix B. 

2. Einstein-Cartan theory 

In the EC theory, the effect of the spin density tensor is locally to induce torsion in the 
structure of space-time. In holonomic coordinates, the torsion tensor Q x fll/ is defined as 
the antisymmetric part of the affine connection T x , 



which vanishes in GR since the connection is assumed to be symmetric in its two lower 
indices. Note that the tilde denotes an EC geometrical object to differentiate it from 
an effective GR object. In the following, Greek indices refer to a holonomic coordinate 
basis, while Latin indices refer to an arbitrary non-holonomic orthonormal basis. 

In order to find a proper description of a Weyssenhoff fluid, we first have to 
determine the EC field equations. The gauge group associated with the EC theory is 
the Poincare group [3] . This is easy to understand as the asymmetry of the connection 
requires an affine generalisation of the Lorentz group which is precisely the Poincare 
group. In the Poincare gauge theory of gravity, the gravitational field is described by the 




(1) 
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tetrad field e^ a and the local spin connection uu ab ' . The spin connection is antisymmetric 
in its Latin indices, u ab = —dJ ba if V \g^ u = which we assume throughout, and 
the inverse of the tetrad is given by e M a , such that e M a e M fe = 5 b and e^ a e v a = 
The geometrical structure of U± — i.e. the metric g^ u and the EC connection T x — 
is completely determined by the tetrad (translational field) and the spin connection 
(rotational field) according to, 

Qtxu = e^e^Vab , (2) 



r p V — e a u) foC^ + e a d u e^ . (3) 

Using the gauge relations (j2J) and ([3]), the torsion tensor (CQ) can be rewritten in terms 
of the translational and rotational fields, 

Q\u = ex a t A [H = -i («9 M e/ - ^e/ + u\^ej> - Q\ v e b ) . (4) 

The metric and the connection are assumed to be compatible, which means that the 
nonmetricity vanishes and implies that the EC connection T x fll/ can be decomposed in 
terms of the Levita-Civita (torsion free) connection Y x and the contortion tensor K x 
as, 



r A ^ — r A ^ - k x ^ u , (5) 



where, 



TV = y Xa ( d ^9au + d u g^ a - d a g^ v ) , (6) 

K ixv = ~Q ixv ~ Qfiv ~ Qvfj, ■ (7) 

The curvature is described by the Riemann-Cartan tensor and its contractions, i.e. the 
Ricci-Cartan tensor and the Ricci-Cartan scalar, 

#V = d^\u ~ du^ bll + ^ c ^% - ^ c hl fi\ v , (8) 

R[a> = R = 6 aCfi R buv 1 (9) 

K = R°\ v = e\e\rf b R a h(JV . (10) 

The field equations of the EC theory are derived from the action S defined on a space- 
time manifold M. as, 



s = f d 4 x \— (n - 2a) + c m 

Jm 12k v ' 



where k = 87rG/c 4 , e = det(e M a ), A is the cosmological constant and C m = 
C m {.^^ a , & ab ^ 4>m) is the Lagrangian density of the matter fields 4> m . Varying the action 
(fTTj) independently for e^ a and Cj ab ^ the field equations are respectively found to be, 

R» a -y a K + e> 1 a A = Kf» a , (12) 

Q\ h + 2e\Q b] = KS» ah , (13) 

where Q a = Q^ afJ- is the torsion trace, and the material sources of the gravitational field 
are respectively the energy-momentum and the spin density tensors defined as, 

f ^ IJCm (14) 
= -^f . (15) 
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These source terms are the functional tensors of the EC classical field theory obtained by 
variation of the action S. They should not be confused with the corresponding canonical 
tensors derived from Noether's theorem since these two kinds of tensors may differ in 
an EC framework. The translational field equation (|T2|) can be recast in terms of purely 
holonomic coordinates and decomposed into symmetric and anti-symmetric parts, 

R{pu) - + 5W A = K T{^v) , (16) 

-R[H = kT^] . (17) 

3. Weyssenhoff fluid description 

The Weyssenhoff fluid is a continuous macroscopic medium which is characterized on 
microscopic scales by the spin of the matter fields. The spin density of matter is 
described by an antisymmetric tensor, 

Sfj,u = ~S U ^ , (18) 

and has been postulated by Obukhov and Korotky [5] to be related to the source of 
torsion according to, 

S\ u = u x S^ , (19) 

where u x is the 4- velocity of the fluid element. The Frenkel condition requires the 
intrinsic spin of a matter field to be spacelike in the rest frame of the fluid, 

S^u" = . (20) 

This condition arises naturally from a rigorous variation of the matter Lagrangian C m 
as shown in [5]. 

The Frenkel condition implies that the torsion trace vanishes, and hence the 
rotational field equations (1131) reduce to an algebraic coupling between spin and torsion 
according to, 

Q\ u = ku x S, u . (21) 

Thus, the torsion contributions to the EC field equations are entirely described in 
terms of the spin density. It is useful to introduce a spin-density scalar S defined as, 

S 2 = \S^S» U > . (22) 

Using the Frenkel condition, Obukhov and Korotky showed [5] that the symmetric 
part of the EC field equations for a perfect fluid with spin (fT6|) can be recast in 
terms of effective GR Einstein field equations with additional spin terms, whereas the 
antisymmetric part f|T7|) simply becomes a GR spin field equation. 

The former are found to be, 

Rfj,u ~ \9y,vR- — K T^ V , (23) 
where the effective stress energy momentum tensor of the fluid is given by, 

T% = {p s + PsWu v - Ps g^ - 2 (g px + u p u x ) V p [u {fl S u)x ] , (24) 
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with effective energy density and pressure of the form, 

p s = P — kS 2 + k~ A , 

p s = p — kS 2 — k~ A , 
satisfying the physical equation of state, 



(25) 



p = wp , 



(26) 



where w is the equation of state parameter. 
The spin field equation is given by, 



(27) 



4. Weyssenhoff fluid dynamics using a 1+3 covariant approach 

We will now use the 1 + 3 covariant approach, outlined for convenience in Appendix A, to 
describe accurately the dynamics of a Weyssenhoff fluid in GR on all scales and in a non- 
perturbative way. Once the dynamical evolution is entirely determined, a perturbation 
analysis can be performed for any given class of models according to their symmetries. In 
a cosmological context, we would require the cosmological fluid to be highly symmetric 
on large scales but allow for generic inhomogeneities on small scales. This is necessary 
to provide an accurate enough description of the observable universe accounting for its 
homogeneity and isotropy on large scales as well as for all the complicated structures it 
contains on small scales. 

In GR, the Weyssenhoff fluid dynamics is actually a generalisation of the dynamics 
of a perfect fluid, where the effective energy density p s and pressure p s contain a 
spin density squared S 2 correction term, and the stress energy momentum tensor T^ u 
incorporates an additional spin divergence term. The new contribution to the effective 
dynamics comes from the spin field equation ( 1271) . Thus, the dynamics of a perfect fluid 
is recovered for a vanishing spin density. 

The dynamical model of a perfect fluid with spin is fully determined by its matter 
content — including the spin properties of the particles — and its curvature. The matter 
content of the Weyssenhoff fluid is described by the effective stress-energy momentum 
tensor ( [241 . Using the 1 + 3 formalism, it can be recast as, 



The physical interpretation of the Weyssenhoff fluid now becomes more transparent. The 
terms containing the effective energy density p s and pressure p s represent the behaviour 
of an effective perfect fluid, where p s and p s account for the spin contributions. The other 
terms describe how the peculiar acceleration of the fluid and the fluid anisotropies 
— described by the rate-of-shear a pv and the vorticity uj^ v respectively — couple to the 
spin density and contribute to the effective energy density of the fluid. 




(28) 



2u^D S u )\ + Au^a S v )\-2a^ S v )\ + 2ojf JJ , S v )\ . 



Weyssenhoff fluid dynamics in general relativity using a 1+3 covariant approach 



7 



All the information related to the curvature is encoded in the Riemann tensor which 
can be decomposed as [8], 

R P » „a = C" vX - 8 p [x R» u] - 5%R" X] - \K8 p [v 5 p x] , (29) 

where C pp vX is the Weyl tensor constructed to be the trace- free part of the Riemann 
tensor. 

By analogy to classical electrodynamics, the Weyl tensor can be split relative to u p 
into an 'electric' and a 'magnetic' part [8] according to, 

= C wua u p u° , (30) 
Hpv = *Cf jL p Pa u p u a = ^t]^ aX C crX vp u p , (31) 

where *C^ vpa is the dual of the Weyl tensor. These parts represent the 'free gravitational 
field', enabling gravitational action at a distance and describing tidal forces and 
gravitational waves. 

The Ricci tensor R pv is simply obtained by substituting the expression (|28|) for the 
effective stress energy momentum tensor T^ u into the Einstein field equations ( J23l) . 

Rfiu = «{| \Ps + 3p s + 8u x S x ) UyUy -\(p s - p s ) 

- 2u {ll D x S u)x + 4u {ll a x S u)x - 2a {il x S u)x + 2u {fl x S v)x } . 

The Riemann tensor R pp vX can be fully split in a 1 + 3 manner according to (1291) 
by using the expression fl32|) for the Ricci tensor R^ u and the decomposition of the Weyl 
tensor C pp vX into its electric E^ u and magnetic H^ u parts. For convenience, the tensor 
is split into three parts: the spinning perfect fluid part (P), the electric part of the Weyl 
tensor (E) and the magnetic part of the Weyl tensor (H). The decomposition yields, 

11 vX — n-P uX i n E vX i* 11 H uX 5 



where 



Rpux = |« {ps + 3 Ps + 12uj x S x ) h [p { X ] u X] - \K Ps h {p [v h" A] 

- 2k (h [p {lJ - u [p u [v ) [-u p] D°S x]a - u x] D a S^ + 2u^a°S x]a + 2u x] a a S^ 
- a^S x]a - a x]a S^ + u^S x]a + uj x]a S^} , 
K»x =C p / uX = AuK [u E^ x] -Ah [p [u E p] x] , 
K\x =C p h \ x = 2 V pp °u [u H x]a + 2 VvXa u^H^ . 

Note that for a vanishing spin density (i.e. in absence of torsion), we recover 
Ellis and van Elst's results p2] after reexpressing the physical variables in terms of 
the opposite signature (—,+,+,-1-). This is also the case for every propagation and 
constraint equation describing the dynamics of the Weyssenhoff fluid because these 
expressions are projections of effective GR identities which are based on the Riemann 
tensor and its contractions. 

In general, there are four sets of dynamical equations for a perfect fluid with spin. 
These sets are derived respectively from the Ricci identities, the Bianchi identities, once- 
and twice-contracted, and the spin field equation. We now discuss each set in turn. 
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4-1. Ricci identities 

The first set of dynamical equations arises from the Ricci identities for the vector field 
u p defining the worldline of every matter field, i.e., 

2V [M V^ P = R [Hp x ux . (34) 

To extract the physical information stored in the Ricci identities, the latter have to be 
projected along the worldlines u p and on the orthogonal spatial hyper surf aces h+ v . The 
non- vanishing projections yield the propagation equations and the constraint equations 
respectively, 

u a h%h\ (2V [a Vp ]U , - R [aP] °u a ) = , (35) 

V pX v h a X x h\ (2 V [o V p] u, - R [af}]j °u a ) = , (36) 

where the latter have been expressed in terms of rank-2 tensors by duality (r] px v ) without 
loss of information. 

The Ricci identities can be further split by separating the propagation and 
constraint equations into their trace part (T), symmetric trace-free part (STF) and 
antisymmetric trace-free part (ATF). The sets of equations are explicitly determined 
by the kinematics of the 1 + 3 covariant formalism (1A.14|) and by substituting the 
Riemann tensor decomposition ( 1331) into the projections yielding the propagation (|3~5l) 
and constraint fl36l) equations respectively before splitting them into parts. 

The propagation equations are found to be as follows. 

• The Raychaudhuri equation (T), 

= -|6 2 + D x a x + 2{u 2 -a 2 - a 2 ) - § {p„ + 3p s + 8u x S x ) , (37) 

which is the basic dynamical equation of a perfect fluid with spin in this system. 
The last term on the RHS describes how the interaction between the spin density 
and the vorticity density affects the large scale dynamics. The physical meaning of 
this term is clear: the energy required to align the spin with the vorticity will act 
like a brake on the expansion, leading to the presence of this damping term in the 
Raychaudhuri equation. 

• The vorticity propagation equation (ATF), 

^<M> = "I + I ( CUrl + , (38) 

which shows how vorticity conservation follows for a perfect fluid. Note that there 
is no spin contribution, which means that torsion does not explicitly affect the 
vorticity evolution, although the effect of spin on the other dynamical variables 
must be taken into account. 

• The shear propagation equation (STF), 
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which shows how the tidal gravitational field and the spin density S^ v induce 
shear. The coupling between the spin density and the shear density contributes to 
the fluid anisotropies by increasing the rate of shear whereas the coupling between 
the spin density and the vorticity density has the opposite effect. 



The constraint equations are given by the following relations. 

• The vorticity divergence constraint (T), 

D x u x = -a x uj x . (40) 

This constraint simply expresses the fact that, in presence of a peculiar acceleration 
induced by a non-gravitational force due to the fluid dynamics, the spatial variation 
of vorticity is proportional to the vorticity. 

• The shear and spin divergence constraint (ATF), 

D x (a/ + uo x + kS x ) - §1^0 = 2a A (uo^ x + kS X ) . (41) 

Using the vorticity constraint (HOI) , the shear and spin density constraint (jUJ) can 
be recast as, 

D x (a/ + kS x ) - |D M 9 = a x (3a;/ + 2kS x ) . (42) 

This expression relates the spatial variation of physical quantities, such as the spin 
density, the rate of shear and the expansion rate on the LHS, to the coupling 
between the acceleration due to the fluid dynamics and the fluid anisotropies on 
the RHS. 

• The magnetic constraint (STF), 

= -DfrU„) + 2a^u v) + (curl o)^ . (43) 

Using the vorticity constraint (|4"0|) . the magnetic constraint (j4"3l reduces to, 

Hfj, u = 3a^u u) + (curl a)^ . (44) 

This constraint shows that the magnetic part of the Weyl tensor is induced by 
the curl of the shear and the coupling between the acceleration due to the fluid 
dynamics and the vorticity. 

4-2. Once- contracted Bianchi identities 

The second and third set of dynamical equations are contained in the Bianchi identities. 
The Riemann tensor satisfies the Bianchi identities as follows, 

V^R Xv \ p = . (45) 

By substituting the splitting ( 1291) of the Riemann tensor R Xu ^ p and the effective 
Einstein field equations fl23|) into the Bianchi identities (j4"5l) and contracting two indices 
(a and p), the once-contracted Bianchi identities are found to be, 

V p C% p + V [X R U \ + [x V u] 7Z = . (46) 
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In a similar manner to the Ricci identities, the information stored in the once-contracted 
Bianchi identities has to be projected along the worldlines and on the orthogonal 
hypersurfaces h^ u . The projections yield respectively two propagation and two constraint 
equations, 

h a{ ,h\ )Up (Vc a % + V [a R?\ + [Q V^) = , (47) 

Vx^K^h^ + V^ a R% + i<J 7 [a V%) = , (48) 

K,h\ (VC^ + V^ 7 + i<5 7 ' a V%) = , (49) 

T] X ^v?h\h^ (v p C a % p + V [a R p \ + |5 7 ' a V'7l) = . (50) 

The sets of equations are explicitly determined by substituting the expression for 
the Weyl tensor splitting fl33|) and the Ricci tensor fl32|) into the projections of the 
once-contracted Bianchi identities ( 1471) — ( 1501) . 

The propagation equations are found to be as follows. 

• The electric propagation equation, 

E(llv) = - @ Efa, + (curl H) pu - § (p s + p s ) 



+ 3a {p x E u)x + u {tl x E u)x - 2ri pX ^H v) x a p + k (5^)^ , 



(51) 



where 



— §cr Ap (o r < At A 5'i/) p — ^-V A >SV} P ) 

+ I (D {ll - 2a (M ) (£>% A - 2a% A ) . 

This equation is similar in form to Maxwell's electric propagation equation in an 
expanding universe. The {S^),, term on the RHS of relation (15TI) describes how 
the coupling between the spin density and the fluid anisotropies contributes to the 
gravitational tidal field E^. 

The magnetic propagation equation, 

H {IMU) =-QH, u - (curl E)^ 

A A A J 

+ 3(7^ H v )x — ^(n Hv)\ + 2rj p x{fj,E u ) a p + k(S^)^ , 

where 

OStfV) = f^U^K* - " pX )S„)x + S p \a v)x - u v) x)} 

- (a v f - u v) °){D x S pX -2a x S px ) 

- aj ap (D x S v)x - 2a x S u)x )} . 

This expression is analogous to Maxwell's magnetic propagation equation in an 
expanding universe. The (Sh)u v \ term on the RHS of this relation (152!) describe 
how the coupling between the spin density and the fluid anisotropies contributes to 
the gravitational tidal field H pv . 
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In a similar manner to that in which Maxwell's equations describe electrodynamics 
in an expanding universe, the coupling between the electric (jSTj) and magnetic fl52l) 
propagation equations gives rise to gravitational waves damped by the expansion of the 
universe. 

The constraint equations are given by the following relations. 

• The electric constraint equation, 

D x E pX = § D pPs - ?>u x H pX - v^xo- u p H x ? + k (S divB ) M , (53) 

where 

(S divE ) M = D fl [((T Xp - ^)S Xp } - D x [{a { ; - co { ;)S X)p ] 

- |6 (d x S, x - 2a x S pX ) + \a w (D X S<> X - 2a x S pX ) . 

This constraint is a vector analogue of the Newtonian Poisson equation. It is 
similar in form to Maxwell's electric divergence equation. For this gravitational 
field equation, the source is not the electric charge density but the energy density. 
The (SdivE^ term on the RHS of expression (153]) describes how the coupling between 
the spin density and the fluid anisotropies acts like an effective electric divergence 
source. 

• The magnetic constraint equation, 

D x H pX = k{ Ps + p s ) uj p + 3u x E pX + rj^xa^P + « (S div ^ , (54) 

where 

(S divH ) p = \v^ P D v (D X SP X - 2a x S" x ) . 

This constraint is analogous to Maxwell's magnetic divergence equation. Unlike for 
Maxwell's equation, this gravitational field equation has a source term which is the 
fluid vorticity. The (S , di v if)„ term on the RHS of expression flMl) describes how the 
coupling between the spin density and the fluid anisotropies acts like an effective 
magnetic divergence source. 



4-3. Twice- contracted Bianchi identities 

The third set of equations is given by the twice-contracted Bianchi identities which 
represent the conservation of the effective stress energy momentum tensor. They are 
obtained by performing a second contraction (fi = v) on the once-contracted Bianchi 
identities fT4"6T) . 

v" (iv + \g pv n) = kvt; = o . (55) 

There are only two possible projections to extract the information stored in the 
twice-contracted Bianchi identities, 

u^T% = , (56) 
hfVTl, = . (57) 
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The propagation and constraint equations are explicitly determined by substituting 
the reduced expression for the stress-energy momentum tensor (1281) into the two 
projections of the twice-contracted Bianchi identities ( 1561) and ( 1571) respectively. 

The propagation equation is found to be as follows. 

• The effective energy conservation equation, 

p s = -6 {ps+Ps) . (58) 

Note that for a vanishing spin density this relation reduces to the well-known energy 
conservation equation determining the evolution of the physical energy density p 
and pressure p. 

The constraint equation is given by the following relation. 

• The momentum conservation equation, 

D^Ps = (p s + Ps) ctf, + (Sp)^ , (59) 



(S p )^ = — 2 (D u — a") (a^S^x - oj { „ x S v)x ) - (D x S II x - 2a x S flX } 
- |6 ( D X S„, - 2a x S„,) - (a„ v - uj„ v ) ( D X S^ - 2a x S„ 



where 



go o^x - ^« J^x J \y(jt "'n ) \v o vX - zu D v Xj 

The term (S p ) describes how the coupling between the spin density and the fluid 
anisotropies contributes to the total angular momentum. 

4-4- Spin dynamics 

The last dynamical equation for the evolution of the Weyssenhoff fluid is the spin field 
equation (T27j) . To extract the spin propagation equation, the field equation has to be 
twice projected on the hypersurface orthogonal to the worldline. By duality, we can 
write it in terms of the spin density pseudovector without loss of information (1A.19j) . 
and we obtain: 

• The spin propagation equation, 

£ w = -es M . (60) 

This expression ( 160|) can be recast in terms of the spin-density scalar S 2 ( 1221) defined 

as, 

S 2 = -S^S" . (61) 

It is then simply given by, 

S = -QS. (62) 

This relation shows that the evolution of the spin density is the same on all scales 
because it is entirely determined by the volume rate of expansion of the fluid. 
For consistency, note that this expression implies that the spin density is inversely 
proportional to the volume of the fluid. 
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The effective energy conservation equation (1581) can now be recast in terms of the 
true (i.e. not effective) energy density and pressure of the fluid by substituting the 
spin propagation equation ( 1621) . 

p = -Q (p+p) . (63) 

The effective energy density p s and pressure p s contain spin density squared S 2 
correction terms (f25l) . Thus, the spin propagation equation (1621) and the energy 
conservation equation (163]) imply that the spin density will rule entirely the dynamics 
of the fluid at early times (kS 2 ^> p,p), whereas, at late times, the spin contribution 
can safely be neglected (nS 2 <C p,p)- 

In a cosmological context, the spin dominated era might lead to an infiationnary 
behaviour. This promising prospect will be analysed in detail in further work. Given 
that the matter dominated era is not affected by the spin contribution, the cosmological 
model thus reduces to the dynamical behaviour of a perfect fluid in GR. Hence, the spin 
density contribution from the Weyssenhoff fluid is expected to affect significantly the 
early time evolution of the fluid leaving the late time dynamics unchanged. Therefore, 
it is not currently promising as a candidate to describe dark energy. 



5. Consistency of the dynamics for an irrotational Weyssenhoff fluid with 
no peculiar acceleration 

The consistency of the propagation and constraint equations can be verified by evolving 
the constraints. This is a tedious but straightforward task. To make the problem 
tractable, we chose to restrict our attention to the class of models for which the fluid 
dynamics is described by an irrotational flow (i.e. = 0) with no peculiar acceleration 
(i.e. a M = 0). This ensures a hypersurface-orthogonal flow and the existence of a globally 
defined cosmic time. If the flow is initially irrotational, it will remain so at later times 

For each space-time slicing, we can now define the curvature tensors entirely in 
terms of the spatial hypersurface orthogonal to the worldline. For this purpose, let us 
define a vector v A , which is orthogonal to the worldline, and an expansion tensor Q Ml/ 
according to 

v\ x = o , = |e v + <v • (64) 

The Ricci identities on the 3-space orthogonal to the worldline can be defined as 

2D [fJi D v] v p = *R fML , p x v x , (65) 

where the 3-space Riemann tensor *R^ P \ is related to the Riemann tensor R pup \ by 

*R Pl iu\ = h a p h? ^h? u h s xRap-fS + G pi ,G M A — QpxQ^u ■ (66) 

The 3-space Ricci tensor and scalar can be obtained by contracting the 3-space Riemann 

tensor with the induced 3-space metric h pu , 

*R pu = h px * Rp^xu , (67) 
*TZ = h^hP x *R P p Xu . (68) 
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Using (lA.14p . (1551) . ([3HD and (16"B"j) . these 3-space curvature quantities can be recast 
respectively as, 

*R P \ X = - \KhP [v h» x]Pt - + 2Q<> [V QK X] 

(69) 

- iKh\ (e x]a s^ + e^s A](7 ) , 

= o-( Mi/ ) + 9 - na {p x S u)x - \h pv [2k Ps - §6 2 + 2a 2 ) , (70) 



*ft = |e 2 _ 2k Ps - 2ct 2 , (71) 

where the last relation is the generalised Friedmann equation expressed in terms of the 
spatial curvature *TZ. 

5.1. Evolution of the constraints 

To determine the time evolution of the constraint equations, we shall follow Maartens' 
approach [H] and generalise his results to the include the presence of spin. For an 
irrotational Weyssenhoff fluid in absence of any peculiar acceleration, the propagation 
equations ([37]), (J32D, fl5H), (TO and flU, denoted by V A = where A = 0, ... ,5 , 
reduce to 

V\ = S M + QSp , (72) 
V 1 = p + e{p + p) , (73) 
= 6 + |9 2 + 2a 2 + f (p s + 3p s ) , (74) 
^% = + 1 0O V + °"< m A(T ^a + ^ - S»)X , (75) 
P% =£ </w> + 6 E^ - (curl H)^ + f (p s + Ps) <V - S<t</E v) x ^ 

+ « (o^S^a) + f ©O^S^A + f 0-\ p a^ X S u ) P — t;D^ l D X S v )\ , 



7>% =#< H +6ff F + (curl £)^ - 3a ( /^ >A 

- %v« P (,D a (<J pX S v)x + a v)x S px ) + f v^ufD^ , 



(77) 



and the constraint equations (J4"2l) . (jSJ), fl53|) . fl54"|) . and fl59|) . denoted by C A = where 
A = 0, . . . , 4 , become 

C° M =D,p s + 2/J A (<%) - a x D^ x + a x D p S Xp 

- S X D X & - \S x D?S Xp + ^»xS%H x ? , 
C\ = D x a pX + nD x S pX - \D p Q , (79) 
C\ v = (curl a)^ - H pv , (80) 
C\ =D X E pX - %D pPs + r\ ilvX o~ v p H xp 

+ kD x (a { /S x)p ) + l&D x S pX - 1a pp D x S px , 
C\ = D x H pX - i] puX a u p E xp - 1 VpuX D»D p S xp . (82) 

The evolution of the constraints C A along the worldlines u p leads to a system 
of equations C A = J rA (C B ), where T A do not contain time derivatives, since 
these are eliminated via the propagation equations V A and suitable identities. The 
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covariant analysis of propagation and constraint equations involves frequent use of a 
number of algebraic and differential identities governing the kinematical and dynamical 
quantities. In particular, one requires commutation rules for spatial and time derivatives. 
The necessary identities are collected for convenience in Appendix C. After lengthy 
calculations the explicit time evolution of the constraints (1791) . (IBUl) . and ( 152]) is found 
to be, 

(<*) i = -eC\ - 2 V /°a a x C\ p - C\ + kC\ , (83) 

(c%) ± = -ec% + ^v-)P cl A' ( 84 ) 
(c\) ± = - lec\ + Wc\ + \h*c\ 

+ r,^H a x C 2 pX - |curl C\ . 

The constraints are preserved under evolution as now briefly explain. Suppose that 
the constraints are satisfied on an initial spatial hypersurface {t = t }, i.e. C \t = 0, 
where t is the proper time along the worldlines. Since C A = is a solution for the initial 
data, it then follows from ([53"]) — (155]) that the constraints are satisfied for all time. 

The time evolution of C° p was not explicitly established because the equation of 
state needs to be specified for this endeavour. Neither was the expression for the time 
evolution of C s p explicitly determined due to the overwhelming algebraic complexity of 
that particular computation. However, it is plausible that the dynamics is consistent 
since the three time evolution equations for the constraints ( 1531 . (151]) and (1551) involve 
all the constraint and propagation equations. This is true with the exception of V 1 . As 
we discuss in detail below, V 1 is not involved in the time evolution of (153]) . (151]) and 
(155]) . However, Obukhov and Korotky have shown j5], using the Frenkel condition, that 
any perfect fluid with spin in the EC theory has an energy conservation equation of the 
form V 1 . This is sufficient to show independently the consistency of V 1 . 

The time evolution of C 1 ^, ( 1531 . involves the propagation equations V° p , V 2 , V 3 ^ 
and the constraint equations C° p , C 1 ^, C 2 pu . It has been determined by using the covariant 
identities fjOi]) and fjC18]) . 

The time evolution of C 2 pu , ( 151j) . involves the propagation equations V s pu , V 5 pu 
and the constraint equations C 1 ^, C 2 pv . It has been determined by using the covariant 
identities (1H77]) and (ICTTO . 

The time evolution of C 4 ^, ( 155]) . involves the propagation equations V° p , V 3 pu , V 4 ^, 
V 5 pu and the constraint equations C 1 ^, C 2 pu) C 3 p , C A p . It has been determined by using 
the covariant identities (jCT3]l . (105]) . (jCTBl . (1091) and dCTTO]) . 

The constraint equations are not linearly independent given that they satisfy, 

C\ = -fcurl C\ - D X C% X . (86) 

The consistency of the constraint equations can be explicitly inferred from relation 
( 156]) as explained below [H]. For any given spatial hypersurface, i.e. {t = const}, 
the linear dependence ( 1551 of the constraint equations implies the constraint C 4 p is 
satisfied provided that the constraints C 1 ^ and C 2 pu are also satisfied. Moreover, the time 
evolution of C 1 ^ and C 2 pu) described by ( 1531 and (154]) respectively, depends explicitly on 
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C° M and C 3 ^. Hence, if we take as determining D X S^\, C 1 ^ as denning D^Q, C 2 ^ v as 
establishing H^ u and C 3 ^ as setting D^p s , the constraint equations are consistent with 
each other because C 4 ., then follows. 

The consistency of the constraints for a perfect fluid in GR with non vanishing 
vorticity and peculiar acceleration has been established by van Elst [15]. Thus, having 
shown that the dynamics of an irrotational Weyssenhoff fluid in absence of any peculiar 
acceleration (u = a = 0) is consistent, it is very plausible — although not proven — 
that this will remain the case in the general case when the vorticity and the peculiar 
acceleration are considered. Hence, in that case, to establish explicitly the consistency 
of the constraints for such a fluid, the coherence of the terms involving the coupling 
between the spin density, the vorticity density and the peculiar acceleration would have 
to be shown respectively. This would be a extremely laborious algebraic task, but it is, 
in fact, quite likely to be true since the consistency of two different particular cases has 
already been established. 

6. Comparison with previous results 

A first attempt to study the dynamics of a Weyssenhoff fluid in a 1+3 covariant approach 
was initiated by Palle [11]. The results we find in this paper disagree, however, with the 
majority of the results derived by Palle, as we now briefly explain. 

In a similar way to our own procedure, Palle based his analysis on the effective 
Einstein field equations for a Weyssenhoff fluid obtained by Obukhov and Korotky [5], 
which are outlined in relation {1} of his publication. As explicitly stated in his work, 
Palle projects the EC version of the Ricci identities determined by Hehl [IB], 

2V [M V^ P = R^ux + 2Q\ u V x u p , (87) 

which are given in relation {4} of his paper to find the corresponding propagation and 
constraint equations. This stands in direct contradiction with the fact that the 1 + 3 
covariant approach used is based on effective GR field equations. 

Moreover, in Palle's work, there is no mention of the antisymmetric part of the 
EC field equations which lead to the spin field equation. It seems unfeasible to provide 
an accurate description of a cosmological fluid with spin without describing the spin 
dynamics. 

Furthermore, Palle chose to neglect the contributions due to the electric and 
magnetic part of the Weyl tensor but did not provide any explanation for this. Indeed, 
the relation {7} he obtained for the shear propagation equation has no tidal gravitational 
field E^ v contribution, and there is no magnetic constraint equation. To describe the late 
time cosmological evolution, it seems indeed reasonable to neglect the contributions due 
to the primordial free propagating gravitational fields which have been damped by the 
cosmological expansion. However, these fields do significantly affect the early dynamics 
and have to be taken into account in a general description of cosmological models. 

Finally, Palle does not determine the cosmological relations derived from the 
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Bianchi identities. Again, these would be very useful to understand the dynamics of the 
early time evolution of cosmological models. 

Palle has recently clarified p2] certain points relating to the approach he followed 
in analysing the cosmological implications of a Weyssenhoff fluid. Several issues still, 
however, remain a concern, as outlined below. 

It is perfectly legitimate to analyse the Weyssenhoff fluid dynamics within an 
EC framework without resorting to an effective GR framework. In such a case, the 
appropriate way to determine the large scale propagation and constraint equations is 
indeed to project the EC Ricci identites (157)1 on the relevant hypersurfaces, which is 
what Palle seems to have done. To achieve this, the EC Ricci identities have to be 
explicitly determined using the effective EC field equations. Our contention is that the 
only effective field equations {1} mentioned in Palle's paper [TT], and used to perform 
the calculations, are the effective GR field equations obtained by Obukhov and Korotky. 
We believe that the GR field equations are incompatible with EC Ricci identities, which 
would thus invalidate the analysis. 

The physical motivation for using GR field equations is that it provides a more 
natural generalisation for the dynamics of a perfect fluid within GR. Although Palle's 
procedure seems inconsistent, we have nevertheless translated his results within a GR 
framework to be able to compare them. To compare explicitly our results with those 
obtained by Palle, note that the torsion scalar Q he uses is related — due to the algebraic 
coupling between spin and torsion — to our definition of the spin density 5* by, 

Q = kS . (88) 

It is now straightforward to see that neither the propagation equations {5 — 7} nor the 
constraint equations {8 — 9} he found agree with our own corresponding results. The 
detailed comparision and analysis can be found in Appendix D. We hope that it might 
clarify this particular issue. 

With regard to the scope of Palle's paper, on large scales, the contribution of the 
tidal forces to dynamics of the Weyssenhoff fluid can indeed be neglected. Hence, the 
Weyl tensor can safely be ignored in his approach, but it was not stated by Palle that 
only the dynamics on large scales were under consideration. It was important to clarify 
this issue because we have considered the dynamical evolution of Weyssenhoff fluid on 
all scales. 

Finally, let us just mention that, as suggested by Palle in [T7], it might indeed be 
more appropriate to consider an N-body simulation to determine the large scale and 
late time dynamics of a Weyssenhoff fluid in a cosmological context, as Palle suggested. 
However, this seems to us to lie outside our study, as we simply considered the evolution 
of such a fluid on all scales and for all times. 

7. Conclusions 

We have used the 1 + 3 covariant approach to determine the dynamics of a Weyssenhoff 
fluid in a non-perturbative and hence completely general manner. This gauge-invariant 
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procedure leads to a consistent set of seven propagation and six constraint equations. 
These give respectively the time and spatial covariant derivative of the set of dynamical 
variables (p, 6, a, u, E, B, S). Compared to the dynamics of a perfect fluid in GR, there 
is one additional propagation equation which is the spin density propagation equation. 
Note that the spin constraint is included in the shear constraint. 

Acknowledgments 

S D B thanks the Isaac Newton Studentship and the Sunburst Fund for their support. 
The authors also thank Anthony Challinor for insightful and stimulating discussions. 

Appendix A. 1+3 covariant formalism 

We will briefly outline the basics of the 1 + 3 covariant formalism introduced by 
Hawking and extended by Ellis to describe the fluid dynamics within GR in a non- 
perturbative way. The aim of this approach is to study the intrinsic dynamics of fluid 
models in a physically transparent manner. This formalism relies on covariantly defined 
variables, which are gauge-invariant by construction, thus simplifying the methodology 
and clarifying the physical interpretation of the models. Furthermore, the form of the 
metric does not need to be explicitly specified and can remain fully general until the 
dynamics is determined. Finally, this approach admits a covariant and gauge-invariant 
linearization that allows linearized calculations to be performed in a direct manner [13J. 

To introduce the 1 + 3 covariant formalism, we follow Ellis and Van Elst's approach 
[T2] using the opposite signature. The approach is based on a 1 + 3 decomposition 
of geometric quantities with respect to a fundamental 4-velocity w M which uniquely 
determines the worldline of every infinitesimal volume element of fluid, 

u" = > "X = 1 > (A.i) 

where r is the proper time measured along the worldlines. In the context of a general 
cosmological model, we require that the 4-velocity be chosen in a physical manner such 
that in the FRW limit the dipole of the cosmic microwave background radiation vanishes. 
This condition is necessary to ensure the gauge-invariance of the approach. 

The 4-velocity u M defines locally two projection tensors in a unique fashion, 

= u,u u U\U\ = U% , U^ = l, U, u u» = u, , (A.2) 

h^ = g^- Ufl u u h\h\ = W v , h^ = 3, V^ = °- ( A -3) 

The first projects parallel to the 4-velocity vector u^, and the second determines 
the (orthogonal) metric properties of the instantaneous rest-spaces of observers moving 
with 4-velocity u^. There is also a volume element for the rest-spaces defined as 

V^p = u x r/ Xll up T)n Up = 77^] , rj^ p u p = , (A.4) 
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where rjx^p is the 4-dimensional volume element (r) Xlxup = ri [Xp u P }, ^0123 = \[\ det g^ u |). 
Note that the contraction of the rest-space volume elements can be expressed in terms 
of the induced metric on these rest-spaces as 

rh^rr" = -3!/^ Q /i>'] 7 = -3!/^>' 7] . (A.5) 

Moreover, we define two projected covariant derivatives which are the time projected 
covariant derivative along the worldline (denoted ' ) and the orthogonally projected 
covariant derivative (denoted _D M ). For any general tensor - v , these are respectively 
defined as 

T"-„... " A VaV"-,,.. . (A.6) 
IhT"- „... = //',//",, ...h\... V,7"-,,... . (A.7) 

Furthermore, the dynamics is determined by projected tensors that are orthogonal to 
on every index. The angle brackets are used to denote respectively orthogonal 
projections of vectors and the orthogonally projected symmetric trace-free part (PSTF) 
of rank-2 tensors according to, 

V M = W v v v , (A.8) 
T {i»>) = (h^ p h u \ - \W V h pa ) T pa . (A.9) 



For convenience, the angle brackets are also used to denote the orthogonal projections 
of covariant time derivatives of tensors along the worldline as follows, 

i}W> = h\v v , (A. 10) 

f<v) = (h^ p h v \ - \h^h pa ) f pa . (A.ll) 

The orthogonal projection of the covariant time derivative of a general tensor T^' u is 
denoted by, 



(•/•"••;,...)•. -//"„...//",,... » A V A v-;,... . (A. 12) 

It is also useful to define the projected covariant curl as, 

(curl T^^^I^T...,/ . (A.13) 

Information relating to the kinematics is contained in the covariant derivative of 
which can be split into irreducible parts, defined by their symmetry properties, 

V^u u = + D^u u = u^ay + ^Qh^u + + , (A. 14) 

where 

• = u u V is the relativistic peculiar acceleration vector, representing the degree 
to which matter moves under forces other than gravity. 

• = D^u^ is the scalar describing the volume rate of expansion of the fluid (with 
H — 1 6 the Hubble parameter). 

• = D^u v ) is the trace-free rate-of-shear tensor describing the rate of distortion 
of the matter flow. 
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M „ = D^u u ] is the anti-symmetric vorticity tensor describing the rotation of matter 
relative to a non-rotating frame. 

These kinematical quantities have the following properties, 

a/ = 0, (A. 15) 

<Vu" = , a vll = <v , ^ = , ( A - 16 ) 

u IJa/ u v = 0, u vli = -(J,,,, , 0^ = 0. (A. 17) 

It is useful to introduce two additional pseudovectors known respectively as the vorticity 
and spin density. These pseudovectors are defined by duality as, 

w A = I^^V => «V = -Vn»\^ x , (A. 18) 

5 a = i^^ ^ S IM , = -r, lil/X S x , (A.19) 

and satisfy 

c^u" = , w^" = , (A.20) 

S'X = 0, S^S^O. (A.21) 

It is also of physical interest to introduce three further scalars which are respectively 
the acceleration, the shear and the vorticity magnitudes defined as, 

\aX > , (A.22) 



2 



a 



o- 2 = \a^ v > , (A.23) 
uj 2 = \uj^ u > . (A.24) 

Appendix B. Transformation of physical quantities under a signature 
change 

The signature convention (+, — , — , — ) we have used throughout this paper is the 
opposite of the one (— , +, +, +) adopted by many authors, such as Ellis and Hawking. 
To facilitate the comparison between results obtained using different conventions, the 
explicit transformations for physical quantities evaluated within the effective field theory 
are given below. 

The metrics, the Levi-Civita tensors and the derivatives transform as, 

9nu > ■> hfiu > hfiu ? flpvxp > VfifXp ■> Viw\ > 'HfivX ■> 
d^d^, V M -> V M , 

The kinematical quantities transform as, 

-> , Up > -Up , a^ 1 -> a M , a M -> -a^ , 

— > — (j^ , a;^ — > — cu^ , cj m — > u;' 1 , — > — . 
The dynamical quantities transform as, 
Rpuxp ^ R[u/Xp i 

CfiuXp * C/ui/Ap i 

T — > T 



Rpv 


' Rpi/ , 




-R, 


Eftv ~ 


-> — E 






Spu 




- 
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It is obvious that rising or lowering indices affects the sign of the transformation for 
any physical quantity since the space-time metric g^ v and the spatial metric h pu change 
sign under such a transformation. 

Appendix C. Covariant identites for an irrotational Weyssenhoff fluid with 
no peculiar acceleration 

It is straightforward to show that the derivatives of the induced metric h pv and the 
Levi-Civita tensor r\p V x vanish, 

Dph^ = , = , (C.l) 

D p v^x = , (rjp U x) ± = . (C.2) 

In this appendix, we consider an irrotational Weyssenhoff fluid (u)p U = 0) with no 
peculiar acceleration (a p = 0). The covariant identities are defined in terms of a scalar 
field /, a vector field V p and three tensor fields, A pu , B pu and Cp V satisfying the following 
properties 

vx = , V/ = A„X = , 

Bpv = , Cpi/ — C^p U ) . 



Using the kinematical decomposition ( 1A.14I) . the identities involving the derivatives 
of the scalar field / are found to be, 

D { pD u] f = , (C.3) 

(DJ) ± = Dj - \QDpf - a, x D x f . (C.4) 



Using the Ricci identities flMj) . the identities involving the derivatives of the vector 
field Vp and tensor field Ap U are given by, 

( D iM>)± =DpV v - \QDpV v - o-p X D x V u + Vux P V x H^ 

- Khp [u V p D x S p]x , 
(D X V X ) X = D x Vx - \QD X V X - o~ p xD x V p - kV»D x S p x , (C.6) 
(D x Ap U ) ± =D x Ap U - \QD x Ap V - a Xp D p Ap V 

+ (V P apA% + Vua P Ap a ) H x p {C.l) 

- K (A\hx { p + Ap a h x[u ) D p S a]p , 
(D x Ap X ) ± =D x Ap X - l&D x Apx - o pX D p A x + VlMTp A\H pX 

— f {A a p + 2Ap a ) D p S ap . 

Using the definition of the curl (1A.13j) and the spatial Ricci identities (|65l).the 
identities involving the derivatives of the symmetric trace-free tensor fields B pu and C pu 
yield, 

V P upC\ (curl£r A = -2C pX D { ^B p] x + \Cp P D x B pX , (C.9) 
D\cm\B)px =\V^ P D V (D x B px ) + ^ vp B\ (|6^ A - E> 
+ \^ P o X po\B^ - l^pO^xS^Bf 



7\V\ 



(CIO) 
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(curlfl)^ = (curL3)^ - §0 (curlS)^ - a^rj^D'B^ 

A Act I ^ J 

Appendix D. Explicit comparision with Palle's results 

To compare our results (BHL) explicitly with the corresponding results obtained by 
Palle, we reexpressed his EC propagation and constraint equations — presumably 
obtained within an EC framework — into a GR framework using the relations given 
in Appendix B. The correspondence between the EC and GR connections is given by, 

rV = r V + K («X + «A A + «^/) , (d.i) 

and necessary to recast the EC covariant derivative V M in terms of its GR counterpart 
V M . To be consistent with Palle's procedure, we only considered the dynamics on large 
scales, hence neglecting the contribution due to the tidal forces (E^ = H^ v = 0). 

The propagation equations are respectively found to be (where we highlight in bold 
face the terms that differ): 

(Palle) : 6 = -§0 2 + D x a x + 2 (u 2 - a 2 - a 2 ) - § (p s + 3p s ) , (D.2) 
(BHL) : = -§0 2 + D x a x + 2 (uj 2 - a 2 - a 2 ) - § [p s + 3p s + 8u x S x ) . (D.3) 



(Palle) : d> w = -§0 + \ (curl a) ^ + (a x + kS x ) lo x , (D.4) 
(BHL) : u w = -|6^ + \ (curl a) M + <j x uj x . (D.5) 



(J(^) = - |0 <V + D kii a v) - a^a v) - <7 (/i V„>A + 
Palle : , „ „ N D-6 

+ 2k* { x S„ )x - \h^ v (0 - 2D x a x + 2a 2 + 2k 2 S 2 ) , 

o^) = - |0 <V + ^a*,} - a^a„) - o^V^a + w<^„> 
(BHL) : /a \ ( U -'J 



The constraint equations respectively yield (where we highlight in bold face the 
terms that differ): 

(Palle) : D x (co x + kS x ) = -a x (u x + kS x ) , (D.8) 
(BHL) : D x lo x = -a x u x . (D.9) 



(Palle) : D x (a x + u x ) - f/^0 = 2a x u x , 

(BHL) : D x (a x + oo x + kS x ) - \Dfi = 2a x (u x + kS x ) 



(D.10) 
(D.11) 
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Note that the shear propagation equation is by definition trace free. This result 
is recovered by BHL but not by Palle. Furthermore, in absence of torsion — i.e. for a 
vanishing spin contribution — the shear evolution equation obtained by Palle does not 
reduce to Hawking and Ellis' result whereas the relation obtained by BHL does. 

We could not rigorously verify Palle's result by evolving the constraints because 
the spin contribution to the Bianchi identities are needed for that purpose as shown in 
Section [51 However, it would be of considerable interest if Palle could emulate BHL 
and demonstrate that his set of equations also reduce to Hawking and Ellis' results in 
absence of torsion, and that the consistency of his equations could be established in the 
absence of voritcity and of any peculiar acceleration. 
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